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II libro della natura 6 scritto in lingua matcmatica. 
Galileo Galilei, [II Saggiatore, 1623]. 

"...At the present time, the intimate connection between causality 
and the analytic continuation is revealed. So, it is not improbable 
to develop of a subtraction procedure even in the most general 
case by the use of analytic continuation techniques." 
O.S. Parasiuk, [[7], p.566, the last paragraph, 1956]. 

This possibility is realized explicitly and efficiently in a body of 
our self-consistent renormalization (SCR). By the self-consistency 
is meant that all formal relations between UV-divergent Feynman 
amplitudes are automatically retained as well between their regu- 
lar values obtained in the framework of the SCR. Self-consistent 
renormalization is efficiently applicable on equal grounds both to 
renormalizable and nonrenormalizable theories. SCR furnishes 
new means for the constructive treatment of new subjects: i) UV- 
divergence problems associated with symmetries, Ward identities, 
and quantum anomalies; ii) new relations between finite bare and 
finite physical parameters of quantum field theories. The aim 
of this report is briefly to review main ideas and properties of 
the SCR and clearly to describe three mutually complementary 
algorithms of the SCR that are presented in the form maximally 
suited for practical applications. 



1. Introduction 

The keystone idea on a purely mathematical genesis of 
the ultraviolet (UV) divergencies of the Feynman ampli- 
tudes (FAs) in quantum field theories is at the heart of 
the Bogoliubov-Parasiuk R-operation [TJ H [3J H G3 El E] . 
Using this idea along with related considerations of 
mathematicians of the 19th and 20th centuries,^ the 
author has developed an universal, high-efficient, and 
self-consistent renormalization (SCR) technique which 
is applicable for any dimension n — 2r n + S n ,5 n — 0,1, 
r n e {0 U N + } of a space-time that is endowed by the 
pseudo-euclidean (p, q) metric g^ v , where p + q = n, and 
for an arbitrary topology of Feynman graphs. 

Algorithmically, the SCR is an efficient realization 
of the Bogoliubov-Parasiuk R-operation as some special 
analytical extension of the UV-divergent FAs in two pa- 
rameters lo g and v G by means of recurrence, compati- 
bility, and differential relations fixing a renormalization 
arbitrariness of the R-operation in some universal way 
based on the mathematical properties of FAs only. The 
parameters lo g and v G are depended on a space-time di- 



x It is appropriate to pointed out here that the first regularization recipe of infinities subtraction for turning divergent inte- 
gral into convergent one had been used in Cauchy's "extraordinary integral" [9] 1101 as well as in d'Adhemar's 12 13 and 
Hadamard's 1 141 1151 1161 I17| "finite part of divergent integral". These recipes are similar but not identical. But in both cases it 
was extendeed the validity of the usual rules of change of variable, integration by parts, and differentiation with respect to the upper 
limit of integration to these new objects. The Cauchy's "extraordinary integral" has been used for an efficient analytic continuation of the 
r(z)-function to some noninteger real values Re 2 < firstly by himself Cauchy 1 101 in 1827, and then in the strips (— n — 1 < Re z < — n) 
by Saalschiitz 1181 1191 in 1887-1888. The term "finite part of divergent integral" was introduced by d'Adhemar in his thesis presented at 
the Sorbonne in December 1903, and defended in April 1904, see [ |20|. p. 477]. Referring to Hadamard's article 114) . d'Adhemar [ |13|. 
p. 371], writes "...Independently of each other, we understood the role of these finite parts...". In d'Adhemar's thesis and articles this 
notion was applied to the construction of solutions of equation for cylindrical waves 1 121 113| , whereas Hadamard use finite parts for the 
solution of the Cauchy problem for second order equations with variable coefficients |14| 1151 1 16 1 and an arbitrary number of independent 
variables 1171 . On the applications of d'Adhemar's and Hadamard's "finite part of divergent integral" in greater detail see Hadamard's 
book 1211 . After 40 years later, when analysing the connections between the intuitive and logical ways of mathematical inventions, 
Hadamard I22| wrote: "...All mathematicians must consider themselves as logics. For example, I have been asked by what kind of 
guessing I thought of the device of the "finite part of divergent integral" , which I have used for the integration of partial differential 
equations. Certainly, considering in itself, it lookes typically like "thinking aside". But, in fact, for a long while my mind refused to 
conceive that idea until positively compelled to. I was led to it step by step as the mathematical reader will easily verify if he takes the 
trouble to consult my researches on the subject, especially my Recherches sur les solution fondamentales et I'int'egration des 'equations 
lin'eaires aux d'eriv'ees partielles, 2nd Memoir, especially p. 121 and so on (Annates Scientifiques de V'Ecole Normale Sup'erieure, 
Vol. XXII, 1905) 1161 . I could not avoid it any more than the prisoner in Poa's tale The Pit and Pendulum could avoid the hole at 
the center of his cell...", see [ 1221 . p. 110, and p. 104 or p. 86 in two identical Russian translations from French edition of 1959]. About 
further developments see M. Riesz [23] [24], F. Bureau [25], R. Courant [26], and S.G. Samko, A. A. Kilbas, and O.I. Marichev [27] , 
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D.A.Slavnov [52] that combinatorics of the R-operation 
is overcomplicated considerably and can be simplified 
essentially. 

• Renormalization arbitrariness of the R-operation 
is fixed in such a way that basic functions {R"F) S j = 
(R"F) S j(uj; M c , A) of renormalized FAs obey the 



recurrence relations as basic functions J 7 . 



J-"sj(a;; M e , A) of convergent or dimensionally regularized 
FAs: 



M, J" s _2j-i - + {uj+ j) F Bj = 0, 



+ (u> + j) {R V T) S 



0. 



(2.1) 



The explicit form of T S j and (R^ J 7 )^- are given be- 
low by Eqs.(3.30)-(3.31). On the self-consistent version 
of the Clifford aspect of the dimensional regularization 
which efficiently overcomes the known difficulties con- 
nected with n-dimensional generalization of the Dirac 
7 5 matrix see [53] IM] ■ 

• The compatibility relations of the first kind: 



(R u F) sj = T sh if u sj := [{y - s)/2] + j < -1, 



mcnsion n, a graph-topological invariant \C\ determining 
a number of independent circuits of a graph G, and two 
FAs characteristics A G and d G . The numbers A G and 
d G determine the maximal degree of polynomials of the 
denominator and the numerator respectively in the in- 
tegrand. As a result, the SCR is efficiently applicable 
on equal grounds both to renormalizable and nonrenor- 
malizable theories that is very important for quantum 
gravity. 

By the self-consistency is meant that all formal rela- 
tions between UV-divergent FAs are automatically re- 
tained as well between their regular values obtained 
in the framework of the SCR. The SCR furnishes new 
means for constructive treatment of new subjects: i) UV- 
divergence problems associated with symmetries, Ward 
identities, reduction identities, and quantum anomalies; 
ii) new relations between finite bare and finite physical 
parameters of quantum field theories. 

The aim of this report is briefly to review main ideas 
and properties of the SCR, see Sect. 2-3, and clearly to 
describe three mutually complementary algorithms of 
the SCR, see Sect. 3-5, which are presented in the form 
maximally suited for practical applications. 

2. The bases and possibilities of the SCR 

2.1 The SCR is an efficient realization of the Bogoliubov- 
Parasiuk R-operation [HI21I31I11E1E1E1IB] which is sup- 
plemented with recurrence, compatibility, and differen- 
tial relations fixing a renormalization arbitrariness of the 
R-operation in some universal way based on mathemati- 
cal properties of Feynman amplitudes (FAs) only. In its 
turn, the BP-approach is rested on an idea that the na- 
ture of UV- divergences is purely mathematical and, per 
se, the R-operation is a constructive form of the Hahn- 
Banach theorem on extensions of linear functionals, see 
for example [28 ] 129 ] 130] . 

2.2 Elaborating this idea the author [3T ] I32 ] I33 ] l34l 
ESI [Ml EHl EE S3 SH H21 S31 SH S3 has obtained the 
high-efficient realization of this renormscheme. In this 
realization: 

• Properties of special functions of the hypergeomet- 
ric type are used essentiallyll 

• Combinatoric is simplified considerably. Our in- 
vestigations confirm the very important assertion by 

2 The connection of particular FAs with the hypcrgeometric functions are well known. See, for example, investigations on analytic 
properties of convergent scalar FAs by using of algebraic topology methods 46 47, 48 , or calculations for needs of phcnomenological 



(R^T) 



8+1, j 



(2.2) 



and the compatibility relations of the second kind: 



-7 r s-2.j-i(^; M e , A) = J" S) j_i(u;; M e , A) = 
= F sj (oj-l;M e ,A), 
{R v F) s -2,j-i{u; M c , A) = {R v F) sj {u - 1 
WF) a , j -i( U ;M e ,A) = (R»- 2 T) sj (oj - 1 



(2.3) 



M e , A), 
Me, A), 



are satisfied automatically. From the first one of 
Eqs.(2.2) it follows that formulae for the regular val- 
ues obtained in the framework of SCR describe uniformly 
both divergent and convergent FAs. 

• The differential relations for !F S j and (RgJ-) S j with 
respect to mass-damping variables /i/ := (mf—iei), I G C, 



■dm 



{R v T)sM 



■oti 



F s j{ijJ - m) 
(R»F) S3 (w-m) 



(2.4) 



physics some classes of FAs by using of differential equation method 49 50 51 . But in our case this connection is established for general 
divergent FAs in any spase-time dimension n and the (p, g) pseudo-euqlidean metric, p + g = n. Apart from, this connection suggests 
some simple method of fixing a renormalization arbitrariness of the Bogoliubov-Parasiuk .R-operation in some universal way based on 
the mathematical properties of FAs only. As a result we obtain the self-consistent renormalization with new valuable properties and 
possibilities. 
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are the same, and the differential relations for ones with 
respect to external momenta k e , e E £, 



,/2] 



= 2 r 



x=0 



^■::.tW 



J 7 ,^ (w — m + x) 



F) s j(0J — m + x) 



(2.5) 



are almost the same. Here d°* = d / d{k ei ) ai , and 

- 4 ei'---°^»( x ) = •^ei'''-°e r ^( x l Q! )^) are special homo- 
geneous polynomials of degree m — 2x in A^J = 
Ag* (a, fc) := l^ ee 5 A eie {a)k^* and of degree x in 
(e'iej-) := As iej (a)fl' <r<CTJ where A ee '(a) are matrix ele- 
ments of the quadratic Kirchhoff form in external mo- 
menta k e , e G £. The polynomials ^ 1 r '.'. IT e ™ {>t\a, k) 
have an algebraic structure of quantities generated by 
the Wick formula, which represents a T-product of m 
boson fields in terms of some set of iV-products of m — 2x 
boson fields with k primitive contractions. Here the 
quantities A g 4 and ( ) play the role of boson fields 
and their contractions respectively. 

• It is essential that T s j and (R^J-) S j as functions of 
two variables M e and A are the homogeneous functions 
of the same degree uj + j. From this it follows that they 
are solutions to the same partial differential equations, 
namely to the Euler equation for homogeneous functions 



[M £ d Me +Ad A -(u + j)] 



0, 



(2.6) 



and to some family of second order equations emerging 
from Eq.(2.6), for example 



,2 

AA 



[M e dl cMt ±(M e ±A)d 2 MtA ±Ad[ 

-(u>+j-i)(d Me ±d A )] r - K,U] 



= 0, 



(2.7) 



that can be again represented as the Euler equation 



[M e 8 ltm +Ad A -(u + j-l)]- 
(d Me ± d A )F a j(u) 



= 0. 



(2- 



So, an important role of the quantities (d Mc ±d A )F S j(ui) 
and (d Me ± d A )(RgJ 7 ) S j(uj) is revealed in our problem. 
After repeating this procedure n + 1 times one obtains 



[M e d„. + AB A - (u + j - n - 1)]- 
(d Me ±d A )T^(u J ~N) 

(d Mm ±d A )(RZ?)"±( U -N) 



(2.9) 



where we define J r ™ ± 



^ (lj - n) := (d Mc ± d A ) N F s:j {uj) and 

n) := (d Me ± d A ) N {R^) s ,j{uj). If n such 

-«±, 



that (w-w+j) < — 1 then both (d Mc +d A )F'^ (u—n) 
and (d Mc + d A ){R v Q F) N s ^{uj - n) = 0. As a result, 
Eq.(2.9) with the p/ws sign is degenerated into the identi- 
cal zero and with the minus sign is reduced to the Euler- 
Poisson-Darboux equation 



d 2 



N 



l)/2/ d 



8M £ dA 



M e - A 



^dAL 



OA J 



{R^) N s f{u:-N) 



= 0. 



(2.10) 



Consistency of solutions to Eqs.(2.9)-(2.10) for differ- 
ent preassigned asymptotic of (R"J-) S j at the vicinity 
of A = leads to relations 



d Mi T sj {u) = -F s j{u - I), 
d A T sj {u) = T sj (u - 1), 
9 Me «J-) SJ -H = -«JF) SJ .( w -i), 

d A (R v F) sj (u>) = (R»- 2 F) sj (u-l), 



(2.11) 



that are also followed from the explicit form of the basic 
functions F S j and {R"F) s j, see Eqs.(3.30)-(3.3I) below. 

2.3 Relations (2.1)-(2.11) manifest a mutual consis- 
tency of asymptotic properties of different terms of FAs 
with respect to external momenta and masses. It is pre- 
cisely these recurrence, compatibility, and differential re- 
lations that are of great importance for investigating of 
symmetries and anomalies problems and for turning of 
developed renormscheme into the self- consistent one. 

Besides, there exist some obvious identities of generic 
nature which are called as the reduction identities 
(RIs) 03 H5 that in another way lead to the recur- 
rence relations (2.1). The simple idea of cancelling of 
equal factors in factorized polynomials in a numerator 
and a denominator of integrands is used in RIs. The 
reduction identities also are of great importance for ap- 
plications as an origin new nontrivial identities. Some 
of them have been used essentially in our investiga- 
tions |32| HBSD S3 El S3 153 ESI 153 ESI- 

2.4 From Eqs.(2.1)-(2.11) and the explicit form of the 
basic functions T S j, {R"F) S j, see Eqs.(3.30)-(3.33) and 
(3.36)-(3.40), implies the following important properties 
of the SCR: 

Algorithmic universality. The SCR is a special an- 
alytic continuation of any FA firstly given by an UV- 
divergent integral. In so doing divergence indices v of 
FA's may be as large as one needs. Hereafter, this con- 
tinuation will be named as the regular (i.e., finite) value 
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of this FA. As a result, the regular values of FA's respect 
certain recurrence, compatibility, and differential prop- 
erties of an universal character and have already been re- 
alized efficiently as convergent integrals. Therefore, the 
calculation of FA's corresponding to rcnormalizable and 
nonrenormalizable theories does not differ for the two 
in the framework of this renormscheme. Actually, the 
problem is reduced to calculations of the characteristic 
numbers, lo, v s j, and X s j determining the basic functions 

(Ro T) a j. 

Separation of problems. The SCR clearly and ef- 
ficiently separate the problem of evaluating regular val- 
ues of UV-divergent quantities of quantum field theories 
from that of relations between bare and physical param- 
eters of these theories, i.e., the SCR realizes in prac- 
tice this the very important potential possibility of the 
Bogoliubov-Parasiuk R-operation. 

Conservation of relations. Any formal relation be- 
tween UV-divergent quantities will be retained also be- 
tween regular values of those if the regular values of all 
quantities involved in this relation are calculated by the 
same renormalization index v (the maximum one since, 
otherwise, we cannot guarantee the finiteness of all terms 
in the relation). So, the SCR is automatically consistent 
with the correspondence principle. As a result, the reg- 
ular values obtained in the framefork of the SCR do sat- 
isfy the vector and axial- vector canonical Ward identities 
(CWIs) simultaneously. 

Extraction of anomalies (quantum corrections). 

In the SCR, owing to the analytic continuation tech- 
nique, quantum anomalies (i.e., quantum corrections 
(QCs) more exactly) are automatically accounted for in 
quantities satisfying the CWIs. More specifically, quan- 
tum anomalies (i.e., QCs) reveal themselves either as an 
oversubtraction effect for a non-chiral case and for the 
chiral limit case (in these cases the Schwinger terms con- 
tributions (STCs) of current commutators are zero) or 
as the nonzero STCs for the chiral case. If necessary, 
the explicit form of quantum anomalies (i.e., QCs) can 
be easily extracted as a difference between two regular 
values of the same UV-divergent quantity calculated for 
proper and improper divergence indices. 

2.5 Algorithmically, the SCR is a union of three ef- 
ficient algorithms of finding: 

i) the convergent a-parametric integral representa- 
tions of renormalized FAs with a compact domain of in- 
tegration of the simplex type and with the self-consistent 
basic functions (R^^sj, s = 0, . . . , d G , j = 0, . . . , [s/2]; 

ii) the homogeneous /c-polynomials V G j{m, a, k), j = 
0, 1, . . . , [s/2], of degree (s — 2j) in external momenta 



k e , e G £, being as a-parametric images of homogeneous 
p-polynomials V G {m,p), s = 0,...,<i G , of degree s in 
internal momenta pi, I £ £; 

iii) the a-parametric functions A(a), A(a,k), 
Y t {a,k), X w (a), 1,1' e C. 

3. Parametric integral representations and 
basic functions of FAs in the SCR 

3.1 From the mathematical point of view any Feynman 
amplitude (FA) associated with an oriented graph G, 

G:=< V,£U£ | e u = 0, ±1, v,- L G V, / G £U£ >, 

in which V is a set of vertices; £ is a set of internal lines; 
£ is a set of external lines; and eu is an incidence matrix 
(i.e., a vertex-line incidence matrix) such that: eu = 
if the line I G £ U £ is nonincident to the vertex Uj G V; 
eu = 1 if the line I G £ U £ is outgoing from the vertex 
Vi G V; eu = — 1 if the line I G £ U £ is incoming to the 
vertex Vi G V, can be always represented by the integral 

(d n P r :=d n Pl ---d n m , d n Pl :=m=idp?, 

I G £, m:= (mi, . . . ,m| £ |), (3.1) 

p:= (pi,...,P|£|), k:= (fci,...,fc| £ |). 

Here V G (m,p) and Q G (m,p) are a numerator and a de- 
nominator polynomials 

P G (m,p) := n„ i£V Pi (m, p)U leC P l K p) = 

= E!> s G Kp), (3-2) 
Q G (m,p) £ :^Y\iec(^-pf) X ', 
file := mf — iei, mi > 0, e; > 0, A; G N+, V/ G £, 

S G (p, k) is a product of vertex 5-functions 

5 G {p,k):^Y\ v ^ v 5 l { Pl k), 

5i(p, k) := S(J2iec e uPl + Eee£ e ^k e ) ; (3.3) 

|^4| is a number of elements of some finite set A; N + is 
the set positive integers; V G (m,p), s — 0,...,d G , are 
s-degrec homogeneous polynomials in internal momenta 
pi, I G £; Pi(m,p) and Pi(m,p) are multiplicative gen- 
crating polynomials of the numerator V G {m,p) that are 
corresponded to the vertex Wi-contribution Vi(m,p,k), 
and to the internal line ^-contribution A;(m,p) e , respec- 
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tively: 

Vi{m,p, k):= Pi(m,p)5i(p, k), 
deg p Pi(m,p) =:di > 0, Mv t G V, 
a / % Pi(m,p) 

AiKp)e:= (^-WF' 

deg p Pi(m,p) =:di > 0, Ml G C. 

The non-degenerate metric form 

diag<^:= (1, 1, -1, -1 ), 
p i 

p + q = n = 2r n + 5 n , 6 n = Q,l, r„G{0UN + }, 



(3.4) 



(3.5) 



is used for each n-dimensional ^-integration in Eq.(3.1) 
3.2 Two characreristics 



2uj c 



d G , 



(n/2)\C\ - A £ 



|C| = |£|-|V| + 1, A £ :=£ ;e£ A,, 



d G := d v +d L = Y,v,ev d i 



(3.6) 



of the integral (3.1) are especially important. Here \C\ is 
a number of independent circuits of the graph G. There 
exist analogous characteristics for all one particle irre- 
ducible (lPI)_subgraphs G C G. If v G > or u° > 
for some 1PI G C G, the integral is UV-divergent and a 
renormalization is needed [55J [3D] . 

While Eqs.(3.1)-(3.3) are identical to the well-known 
representation in terms of vertex-line contributions, 



S G ( P ,k) 



V G (m,p) 
Q G (m,p) e 



II V t (m,p,k)Y[A l (m,p) e , 



iec 



they are more suited for practical calculations. The uni- 
versal decomposition of V G (m 1 p) in terms of s-degree 
homogeneous p-polynomials V G (m,p) is very useful. 

3.3 We shall use of the Fock-Schwinger exponential 
a-representation, see for example, [SSJ|Bni|BT], along with 
the Hepp regularization [30] that introduces parameters 
ri > in the vicinity of ai = 0, Ml G £, 

, 1 2U = Urn r d ^4^e-^-*\ 
pl = {-id/dqi r y^\ l= o, (3.7) 

(pi ■ qi)-= pirqiag Tcr , < n < a t < oo, Ml e C. 

Then the ratio of polynomials V G (m,p)/Q G (m,p) e in 



Eqs.(3.1)-(3.2) can be represented in the form 
Q G {m,p) e n->o ^.,,.1 



c l f-»\ 
(r) 



dv G {a)i Xc 



s=0 



qi=o 
viec 



W qc 

PC 



{pi ■ ateepe) + (pi ■ q c ) 



(3.8) 



lec iec 



In Eq.(3.8) the p L and q L are (\C\ x n)-dimensional actual 
and auxilary internal momenta column-vectors associ- 
ated with the set of internal lines, C, of a graph G; the 
t is the transpose sign, so p T c is the row-vector; the a cc 
is the |£|-dimensional diagonal matrix of a-parameters; 
the Ac is defined in Eq.(3.6). Here, the integration mea- 
sure dv G (a), the integration region ijjf' (r), and the a- 
parametric function M e = M(m, a) e which is the linear 
form in the square of internal masses with ie-damping 
are defined as, 



dv G (a):= Yl 



iec 



dai a 



A,-l 



r(A* 



i?^'(f):= {a/|0<r; < a t < oo, MleC,}, 



OLipu 



fiu := (mf - iei). 



(3.9) 



iec 



Now, substituting of Eq.(3.8) into Eq.(3.1) and inter- 
changing the order of integration in pi and cti, Ml G £, 
we obtain the very useful representation of the regular- 
ized by Hepp integral 7 G (m, fc)£. The integrand of it 
is the (|£| x n)-dimensional pseudo-Euclidean Gaussian- 
like expression but in depended variables pi,Ml G £, 



/ G (m,fc)f:= c c 



dv G (a) 



d G °? 
Y,V G (m,-id/dq c ) / (d n pY8 G {p Cl k £ )i^ 

s=0 J 



e -iM e + 



qi=0 

viec 



(3.10) 



The set of internal lines, C, can be always decompose 
(as a rule in more than one way) into two mutually dis- 
joint subsets, C = M U7V, M. flTV = 0, which determine 
some skeleton tree, i.e., 1-tree subgraph G(V, M. U £), 
with \M\ = \V\ — 1, and corresponding to it co-tree sub- 
graph G(V, Af U £), with |JV| = |£| -|V| + 1= \C\ of the 
graph G. Supports of all Si(p c , fc £ )-functions, Mvi G V, 
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see Eq.(3.3), are defined by Eqs.(3.11) and are equivalent 
to the matrix relations given in Eqs.(3.12) and Sec. 5, 



^imVm + ZjnPn + e j£ k £ = 0, Vj -the basis vertex, 
Pc = ecmPm + e C£ (j)k £ , 

Pm = ZmnPn + e M£ (j)k £ . (3.12) 



Thus, (|.M| x n)-dimensional integration by means of 
Si(Pc, fc £ )-functions, VUj G V/j, (this is equivalent to 
make use of Eqs.(3.12)), gives rise to the following in- 
termediate a-parametric representation, 



I G (m,kf e :=6 a (k e )c c 



dv G (a) 



■Y,'Ps{m,-id/dq c ) / (d n p) 



-iM e + iW q u c t 



91=0 

Vie£ 



(3.13) 



+ (fcj • £ ££ (j|a)fe £ ) + • e C£ (j)k £ ), 
In-= n £Ar (i|a)fc £ + \e T CH q c , 
6 a (k £ ):=5(32 ee£ e(v*) e k e ). 



The explicit forms and some properties of the matrices 
ecu, e ££ (j), Cm(<x), E ££ (j\a), and ±I £Ar (j|a) are given 
in Eqs.(5.1)-(5.4). 

The change of variables p^ by means of the nondc- 
generate linear transformation such that, 



Pa,- = B Ar ^(a)p^ / - B M ^{a)B^ M (a)f N , 
B^ N {a)C MN {a)B MN (a) = l MN , 
B NM {a)B£ N {a) = C£j>(a). 

del B " 



detfl^a) = [detCW(a)]- 1/2 =:A(a) _1 - 
(d n p y = (d n p ) N | det B MM {ot) | ™ | defcff| w 
= (d"p) Af /A(a)"/ 2 , det 5 = (-!)«, 



L/2 



(3.14) 



to reduce Eqs.(3.13)-(3.14) to the form 

dv G (a) 



I G {m,k)::=8 G {k e )c c 



(f) 



A™/ 2 



d G OO 

■Y.V^m^-id/dq,) f (d"py e 

s=0 "L 



i(pJfPAr) j^c 



-ilf, + 



9i=0 

v/e£ 



(3.15) 



== ' Cvi, (a)/ v )+ 
+ (fcj • E ££ {j\a)k £ ) + {q T L ■ e c£ {j)k £ ) = 
= {k T £ ■ A ££ {j\a)k £ ) + 
+ {q T c -Y C£ {j\a)k £ )-\{ql-X cc {a)q c ). 
Taking into account the formula 
f°° dte ±lt2 = n^ 2 e ±l ^ 4 , 

J — OO ' 

which is followed from [55], Ch. 1.5., eqs. (31-32), we find 

J°° d n pie i tf = ■ K n/2 e i( - p ~ q)7T/i = ^"/^(p-™/ 2 ^ 

J™ go (d n pY r e i &M-*"k Xe = (7r n / 2 i p ) |Af| r w . (3.16) 

So, all n-dimcnsional pseudo-Euclidean momenta inte- 
grations in the (p, (j)-metric are performed. Thus, any 
FA (3.1) is led to the a-parametric representation in the 
fully-exponential form, 



I G (m,kf e := (2ir) n 5 G (k £ )a c 



dv G (a) 
Rf(f) A"/ 2 



^V?(m,-id/dq c ), 



-iM,+iWl c 



s=0 



qi=0 



(3.17) 



W**:= (fcj ■ A ££ (j\a)k £ )+ 



+ (q T c ■ Y C£ {j\a)k £ ) - \{q T c ■ X cc (a)q c ), 



a G :=c G (WV) |Af| (27r) 



W\ = \c\, 



where n-dimensional auxilary momenta, qi,l £ C, still 
are used. The explicit form and important properties 
of matrices A ££ (j\a), Y C£ {j\a), and X cc (a) are given in 
Eqs.(5.5)-(5.13). Some properties of them are new. 

Next, the following two operations must be carried 
out: i) to differentiate the exponential function exp{z(gj- 
Y C£ (j\a)k £ ) - i/ 4 (<?c • X cc (a)q c )} by means of the s- 
homogeneous differential polynomials V° {m, —id / dq c ) 
in —id/dqia, I € L, a £ {1, . . . , n}, 1 < s < d G ; ii) to put 
q la =0,\fl e C, Vcr e {l,...,n}, andVs e {0, 1, . . . , d G }. 
Finally, we obtain the important a-parametric represen- 
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tat ion for the general FA (3.1) 



I G (m,kf e := (27T) n 5 G (k £ )b c 



Rf(r) 



dv G (a) 



d G [a/2] 



— u>— j — iMe+iA 



(3.18) 



■^72 2J2J p 5( m ' a »*)* 

s =o j=o 

A = A(a, fc) := (fcj • A ££ (j|a)fc £ ), 

&°:=c°(7r n / 2 i , ') |c| (27r)- n , \C\ = \Af\, 

P G (m, pa, rk) = P -Jt s - 2 J P G (m, a, k). 

Here [s/2] means the largest integer < s/2, i.e., the 
integer part of the s/2; the quadratic Kirhchhoff form 
A(a, k) in external momenta k e , e E £, and the Kirhch- 
hoff determinant A(a) := det C^x (a) are dehned by the 
topological structure of a graph G and are homogeneous 
functions of the first and |C|th degrees in a, respec- 
tively, see Sec. 5 more detail. The quantities V G 2 \m, a, k) 
are homogeneous fc-polynomials in external momenta k e , 
e e £, of the degree s - 2j, j = 0, 1, ... , [s/2]. They 
are a-parametric images of homogeneous polynomials 
V G (m,p). Each monomial of V°j(m, a, k) is aproduct of 
s — 2j linear Kirchhoff forms Yj(a, k) := J2 e e£ Yie( a )k e 
and j line-correlator functions Xw(a), 1,1' £ C, of a 
graph G. Parametric functions Yi{a,k) and Xw{a) are 
homogeneous functions of the 0th and (-l)st degree in 
a, respectively. See Sec. 4 and 5 more detail. 
By introducing the new variables, 



ai=pa[, VleC/j, oij = p(l-Y, 



iec/j a 0> 



lec 



lec/j 



lec 



R ] + l {r) -> R\{r) x E 1 



£1-1 



r > 0, 



(3.19) 



and assuming that r; = r > 0, V/ £ C, we can perform 
the integration over the variable p, < r < p < oo, by 
using [53], see Ch. 6.3., eq.(3) or Ch. 6.5., eq.(29), 



I G (m,k) r e := (2ir) n S G (k £ )b c 



dp G (a) 
A™/ 2 



si^i- 1 



d G [s/2] 
s=0 j=0 



-W-J-lg-tp V e 



J%((j;V e ):=i- u -> J dpp 

r 

= V e u+j r{-u> - j; irK), F e := M e - A, 
cj:= (n/2)|C| - A £ , fo G := c G (W 2 i p ) |C| (2tt)- 



(3.20) 



Here p involved in the 6 G is the number of positive 
squares in a space-time metric g 11 " . The integration 
measure dji G {a) and the integration domain E'^' -1 of 
the simplex type are defined as 



dp G (a):=S(l-Y,ai)l[ 



lec 



lec 



dai a t 1 



E^l -1 — {ai\ai >0,VleC,^2ai = 1}. 



(3.21) 



lec 



In Eq.(3.21) T(a;x) is one in a two incomplete gamma 
functions appearing in the decomposition, T(a) = 
T(a;x) + -f(a;x), see [64], Ch. 9.1., eqs.(l-2), such that 
at Re a > 0, r(a;0) = T(a), 7(a;0) = 0, where T(a) is 
an ordinary gamma function. 

It is useful to remark that we actually have the reg- 
ularization which combines three ones: i) the Hepp reg- 
ularization |30j . (due to the change in the region of in- 
tegration over auxilary variable p); ii) the analytic reg- 
ularization [80], (due to the complexification of the pa- 
rameter X c , half-degree of the denominator polynomial); 
iii) the dimensional regularization [531 154] , (due to the 
complexification of the parameter n, the space-time di- 
mension). Recall, that A £ and n are constituents of u>. 

For convergent FAs the quantities to + j < 0, Vj £ 
{0, 1, ... , [d G /2} and there exists the limit r — > 0. After 
passing to the limit r — > in Eq.(3.20) we obtain, 



I G (m,k) e := (27r) n S G (k £ )b G 



d G [s/2] 



dp G (a) 
A"/2 



• v i ( m > «> M <> A )> 

s=0 j=0 



^(uiM e) A):=r u " j / dpp-"-i- x e- ip ^ M *- A ^ 



M? +i (l-Z e ) u+ iT(-u-j) 
M e ^'f;r(-w-j+i)|, Z e :=A/Af 6 . 



(3.22) 



fc=0 



It is easily verified that basic functions J- S j(uj; M e , A) 
satisfy Eqs.(2.1), (2.3), (2.11), and Eqs.(2.6)-(2.8). 

In the case of divergent FAs for which u) + j > 
at least for one j G {0, 1, ... , [d G /2}, the limit r — >• 
does not exist. In this case, the expressions (3.6r)-(3.6s) 



strictly defined in the region R \ 



R\{r) x El £ l" 



must be made meaningful in a wider region R + 
R\ x E^l" 1 , where R\:= R\{r)\ r = . 
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The Bogoliubov-Parasiuk subtraction procedure ap- 
plied for this purpose replaces I G (m, k) r t by 



(R»I) G (m,k) e = (2ir) n 5 G (k £ )b G 

dv G (a)(R»l) G {m,a,k) e , 



(3.23) 



(R»l) G (m,a,k) e := 

"id 13 

:=l G (m,a,k) e - ^ ~ TOa lG ( m ^ a > rk ^ 



/3=0 



d 



u+l 



dT{l~TY^-^T G {m iai Tk) t 



(3.24) 



where subtraction operations under the integral sign are 
performed by using the Schlomilch integro-differential 
formula, see the 2nd line of Eq.(3.24), for the remain- 
der term of Maclaurin's series. Firstly, this formula was 
applied explicitly to the FAs in the Parasiuk paper [8]. 
Although this expression guarantees a compact represen- 
tation of the subtraction procedure, it is, nevertheless, 
inconvenient for computational purposes, because it in- 
volves additional integration and differentiations in the 
integrand. The expression in the 1st line of Eq.(3.24) is 
all the more inconvenient for these purposes, since every 
term on the right-hand side of one may be assosiated 
with a divergent integral. 

At the same time the algorithm proposed and ap- 
plied in [311 [321 [331 [341 [351 [361 [371 [381 E9] is based on the 
observation, see [51], Ch. 9.2., eqs.(16, 17, 18), that, 



E 



k=0 
k 



X 

A 7 ! 



e x; y(l + v sj ;x), j(a;x) 



j(a;x) 



k=0 

7(0; x) + r(a; x) = 1. 



r(«) 

r(a; x) 
~1>P 



(3.25) 



Now, if we use: the explicit form of the integrand in 
Eq.(3.18); the homogeneous properties for parametric 
functions in k e ,e £ £, see Eqs.(5.8); the 1st line of 
Eq.(3.24); the 1st line of Eq.(3.25); and the following 
useful relation, 



^ P ! Or? 1 ' J I r=0 



,3=0 



fe=0 



fe! 



v aj := [{v - s)/2] + j, 



(3.26) 



we arrive at the multiplicative realization of the sub- 
traction procedure in the integrand of Eq.(3.23) for the 



regular value of general FA (3.1), 



, d G [s/2] 

■V G (m, a, fc)r tJ - J e- iy =7(l + v sj -iA). (3.27) 

The integral in Eq.(3.23) with the integrand (3.27) at 
v > v G is now well-dehned in the domain The 
substitution of (3.27) into the integral (3.23) and the 
change of variables in integration according to Eq.(3.19) 
give rise to the expression 

dp G (a) 



(R»I) G (m,k) e = (2tt)" 6 G (k)b c 



A™/ 2 



ei^i- 1 



d G [s/2] 

■ E E v °i ( TO ' a ' fc ) w ?)>i ( w ; M « A ) . 

s=0 2=0 

{R" F) 8j (w;M t ,A);=i 



■Ul-J 



(3.28) 



dpp- u -*- 1 e- i '* M '- A: > 7(1 + v s j\ipA) 



M? +j q JJ Xsj) , Z e 1 +^ a F 1 (l,A. i ; 2 + !/„•; Z e ), 



l T(2 + u sj ) 

v sj := [(v - s)/2] + j, X SJ :=-u-j + l + 

The integration over p in (3.28) is performed with the 
use of the formula, see [55], Ch. 17.3., eq.(15), 



dxx^ 1 e vx j(v;ax 
a"T(p + v) 



%Fx{l,p+ v\ 1 



(a + v)^ +l T(l + v) " " y ' r " ' " ' " ' a + v 
Re (a + v) > 0, Re v > 0, Re (p + v) > 0. 



3.4 So, using properties of special functions substan- 
tially, the author has obtained [31113113213313313111101 
EU 02 US] high-efficient formulas which realize an ana- 
lytical continuation (in the variables lu g and v G ) of the 
FAs which are represented first in Eqs.(3.1)-(3.3) by UV- 
divergent integrals, and are given finally in Eqs.(3.28)- 
(3.30) as convergent ones. As a result, we have the fol- 
lowing a-parametric integral representation, 



I G (m, k) e 
(R»I) G (m,k) f 



(2ir) n 6 G (k) b c 



dp G (a) 
A™/ 2 



ei^i- 1 



d G [s/2] 

EE p S( m ' a ' fc ) 



F sj (oj;M e ,A) 
(R»F) sj (u;M e ,A) 



(3.29) 



s=0 j=0 

for convergent or dimensionally regularized value 
I G (m, fc) £ , and for regular value (RgI) G (m,,k) e of the 
integral (3.1). The subscript and superscript v on R 
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indicate that (R"I) G (m, k) e is the regular function in the 
vicinity of zero values of external momenta k e , e € £ , and 
is evaluated for an renormalization index v = v G . 

The explicit forms of basic functions F S j(uj;M e ,A) 
and {RgF) s j{oj\M e ,A) are as follows: 

T aj (w; M e , A) := M^{1 - Z e ) T(- uj - j) 

oo k 

= M e u+i 52T(-w-j + k) Z f _ := A/M e , 



k=0 



(R^) sj (oj; M t ,A) := M e w+ ^ r(A„-)/r(2 + !/„•)• 
.Z e 1+ ^ 2 iq(l,A si ; 2 + !/„•; Z e ) 



[(1/ - s)/2] + j = [u] + j + <t s 



'S3 ■ 



\ sj :=-u}-j + l + v sj = 1 - S n 5 lcl /2 + a s , 
[uj] := r n \C\ + 8 n r m - A, w = [uj] + 5 n 5 lcl /2, 
a s := [(5 n 5\c\ +d — s)/2], |C| = 2r lc \ + S [C \, 
v = v G ', ui =uj g , A = A G , d = d G . 

The [{v - s)/2], [(1/ + 1 - s)/2], and [w] in Eqs.(3.28)- 
(3.31) are the integer parts of the (v — s)/2, (y + l — s)/2, 
and uj, respectively. The subscripts (s,j) on F S j and 
(R^J-) S j just mean that these functions are attached 
to the homogeneous fc-polynomials V G j(a,m,k) of the 
degree s — 2j, j — 0, ...,[s/2], in external momenta 
fee, e G £. The latter are a-images of the homo- 
geneous p-polynomials Vf{m,p) of the degree s ap- 
pearing in V G (m,-p), see Eqs.(3.2). The fc-polynomials 
V°j(a, m, k) are constructed by means of a-parametric 
functions Y[(a, k) and X^ (a), I, I' G C. The efficient and 
universal algorithm of building V^Aa, m, k) is presented 
in Sec. 4. The a-parametric functions M e = M(m,a) e 
and A = A(a,k), incoming in Eqs.(3.30) are defined 
in Eqs.(3.9) and (3.18), respectively. The M(m, a) e is 
the linear form in the square of internal masses with ie- 
damping. The functions A(a, k) and Yi(a, k) are known 
as the quadratic and linear Kirchhoff forms in exter- 
nal momenta, k e , e € £■ The function A = A(a) 
is the Kirchhoff determinant, and the Xw (a) are the 
line-correlator functions. The high-efficient and univer- 
sal algorithm of finding a-parametric functions A(a, k), 
Yi(a, k), Xip(a), and A(a) is given in Sec. 5. 

3.5 Investigation of complicated tangle of problems 
associated on the one hand with renormalization meth- 
ods and on the other with conserved and broken symme- 
tries, the Ward identities behavior, the Schwinger terms 



contributions, and quantum anomalies requires of find- 
ing renormalized FAs for different divergence indices. 
For example, amplitudes involved in the Ward identities 
have divergence indices v G and v G + 1. 

Regular values (Rg +1 I) G (m,k) f calculated for the 
renormalization index v G + 1, once again, have form of 
Eq.(3.29) but with another basic functions (i?o +1 J-") s:) -: 

(K+lF) SJ := M e u+j r(A s 1 ,)/r(2 + zA) 



,i ._ 



■Z £ + ^ 2 F 1 (l,A^.;2 + ^ i ; Z e ), 
[(u+l-a)/2]+j=[u]+al+j, 



(3.31) 



AL := - u - j + 1 + vjj = 1 + oi - S n 5 m /2, 



(3.30) al:= [(S n 5 m +d + l-s)/2}. 

In general, {R^ +1 T) sj ^ {R V F) S] , as far as v 1 ^ ^ v sj . 
Difference between them is the important quantity 



(A("+^)JV= (R» +1 F) sj 

Au+x,v) r(A s -, s> 



= -e 



T(2 + u sj ) 



{R v F) sj 



(3.32) 



S3 > 



2r v + 6 V , 



■ = o~ s — °~s =| $v — 8s |i 
2r s + 6 s , v,s£{0UN + }, 



where H + (x) is the Heaviside step function such that 
H + (x) =0, x < 0, H + (x) = 1, x > 0, and S V ,S S := 
v(mod 2), s(mod 2) = 0,1. It is this quantity that per- 
mits to obtain some efficient formulas for calculating of 
the quantum corrections (QCs) (i.e., quantum anoma- 
lies) to the canonical Ward identities (CWIs) of the most 
general kind. For example, to the Ward identities in- 
volving canonically non-conserved vector and (or) axial- 
vector currents for nondegenerate fermion systems (i.e., 
for systems with different fermion masses). Another a 
very useful quantity that is produced by differences 

(A(^lF) SJ := (R.: +2 T) S 



— {R V ^)s-2,j - {Ro J 7 )s-2,j 

= (R v a J") s -2,j - (R u F)sj = 

r(V 



■H + (l + u aj ) 



r(2 



M u+j Z l+u,, 



(3.33) 



is closely related with (A 



(.v+l,v) 



3.6 The expressions given by Eqs.(3.28)-(3.30) have 
two the very important properties. 

Firstly, they describe both divergent and convergent 
FAs in the unified manner. Really, due to properties [62] 
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Ch. 2.8, eqs.(4, 19), i.e., 2 Fi(a,P;a; z) = (1 - z)-^ and 



o I"? 1 , o 2^1 (a, c; z)/r(c) 

c— > 2-/, 2=1,2,... 



(3.34) 



(a)l - l(6), "V-Vi(a + 1-1,6 + 1-1; I;*), 



(« -1)1 

in the case a = 1, 6 = A s j = — u; — j + 
from Eqs.(3.30) and (3.34) it follows 



l-Lc 



(KT) SJ = M e u+j r(- w - j)2Ft{l, -w-j\ i; z e ) 



Taj, if v s 



I, I e N 4 



(3.35) 



i.e., the first relation in Eqs.(2.2). 

Secondly, the basic functions {R"T) S j = 
(R"J-) S j(w;Mt,A) of the self-consistently renormalized 
FAs obey the same recurrence relations as the basic func- 
tions J- S j = T S j(ui; M e , A) of convergent or dimension- 
ally regularized FAs. Really, if the recurrence relation, 
see [62] Ch. 2.8, eq.42, i.e., 

(c - b - 1) 2 Fi (a, b; c; z) + b 2 Fi(a, b + 1; c; z)— 
-(c-l)aFi(a,6;c-l;z) = 0, (3.36) 

between contiguous Gauss hypergeometric functions gFj 
in the case a = 1, b = A S j, c = 2 + v S j, to multiply on 
the quantity M e u+j zl +Va3 T{\ sj ) /T(2 + v sj ) and to use 
of relations i/ 8 _2,j-i = v s j, K-2,j-i(u) = + 1, 



and v 



s,j—~L 



'»3 



1, \ s j-i(u)) = \ s j(u)), then we obtain 



the following recurrence relations 
M t {R v Q F) s -2 d -i - A (RZT) s ,j-i+ 

+ ( w + J )(KJ) SJ = o, 



(3.37) 



between basic functions (RgT) S j = (RgJ 7 ) S j(uj;M e ,A), 
i.e., the second relation in Eqs.(2.1). 

3.7 Transformation formulae (see 62 Ch. 2.1.4, 
eqs.22 and 23) of 2 Fi give rise to the representations: 



(-l)T(X S] )A^ 
T{2 + v s j)M t 



Z, 



■ 2 F X \ i )U + j + l- 2 + v. 
(R»T) aj = (M £ - 



Ze-1 

Z f 



S ' J " z c 
r(A aj ) 



i 



Z 1+V.j 



T{2 + v sj ) 
• aFi (1 + w + j + 1; 2 + v sj ; Z e ) 



(3.38) 



(3.39) 



The equation (3.38) and the behaviour of 2 F\{a,b; c; z) 
in the vicinity z — > 1_ to determine completely the 



asymptotic of basic functions (RgT) S j for A < in the 
vicinity M e — > 0, i.e., the chiral limit: 



Af.-vO (-l)T(A sj - l)A v - 



r(i + ^)M e A ^ 1 ' 



if i/ s j > and A S j 



1 > 0: 



(K^ M ^i^ln(l-^/M e ), 



r(i + v sj ) 



(3.40) 



if v s j > and A SJ - 



1 = 0; 



(R»T) SJ Mc ~°r(-u;-j)(-Ar + i, 



if za^- > and X S j 



1 < or u s j < — 1, 



which is equivalent also to the asymptotic behaviour 
of the basic functions in the case A —> — oo, M e ^ 0. 
From Eqs.(3.30) follows four different series of values for 



A 



■v 



1: 



X sj - 1 = -(5 n 5| C |/2+ 

+ (r«j - r a ) + [(5 n 5 m + 5 d - 8 s )/2], 

X sj - \ = {r d - r s ) - 1/2, 8 n 8 ]C \ = 1&5 S > 5 d 
= (u - r s ) + 1/2, S n S ]cl = 1 & 5 d > S s 
= {rd-r s ), <5„<5| C | = k S d > 5 S 

= {r<i - r 8 ) - 1, S n S lcl = & <5 S > <5 d 

d = 2r d + S d , s = 2r s + S s , S n , S lcu S d , S s = 0, 1. 



(3.41) 



(3.42) 



It is evident that Eq.(3.39) to present a multiplicative 
realization of the subtraction procedure explicitly, 



{R»T) 



T. 



S3 



(S^)s3 



• aFi(l + u + j + 1; 2 + t/, i; Z e ). 



(3.43) 



4. Homogeneous k-polynomials V°j (m, a, k) 
of a-parametric representation of FAs 

4.1 From Eq.(3.7) it is evident that the basic func- 
tions {R^!F) S j and the homogeneous fc-polynomials 
"Pgj(m,a,k) in external momenta fc e , e € £, of degree 
s — 2j, j = 0, 1, . . . , [s/2], are two closely related in pairs 
important universal ingrediants of the SCR representa- 
tion of FAs. The latter are a-images of the homogeneous 
p-polynomials Vf{m,p) in internal momenta pi, I G C, 
of degree s, s = 0, 1, . . . , d G , appearing in the numerator 
polynomial V G (m,p) 7 see Eqs.(3.1)-(3.2). 
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Each monomial of V°j (m, a, k) is a product of s — 2j 
linear Kirchhoff forms Yi(a,k) := J2 e e£ ^e( a )^e an d j 
line-correlator functions Xw(a), 1,1' G C, of a graph G. 
The efficient algorithm of finding these expressions from 
initial homogeneous p-polynomials Vf(m,p) in internal 
momenta pi, I € £, of degree s = 0,1,..., c? G , has been 
elaborated in (3TJ [3l [33l |34]. It resembles Wick rela- 
tions between time-ordered and normal products of bo- 
son fields in quantum field theory. The main steps of 
this algorithm are as following. 

• The polynomials Vf (m,a,k) are determined as 

P°(m, a, k) := Vf (m,p)\ Pl=Yl ( a ,k), 

j = 0, s = 0,l,...,d G , (4.1) 

i.e., by the straightforward substitution pi — > Yi(a,k), 
VZ G C, in the polynomials Vf(m,p). 

• The polynomials Vfj(m, a, k), j = l,...,[s/2], 
have the algebraic structure of quantities generated by 
the Wick formula which represents a T-product of s bo- 
son fields in terms of some set of iV-products of s — 2j 
boson fields with j primitive contractions. In this case 
the linear Kirchhoff forms Yf{a,k) and their primitive 
correlators 



Wick type rule yields 

[a/2] 
3=0 

d£{l>-V2i) d{1) d(s} 

V^f^^,k)~^Y^(a,k) 

j 

•t[{Xw j W (ajff^w'-w), (4.4) 

d(t>) ci(c) 

where the summation in the second Eq.(4.4) is extended 
over all partitions d of 4*\ • • • , £s^) according to the 
Young scheme (1 8 ~ 2 ^). Then the image of homoge- 
neous p-polynomials Vg(m,p) given by Eq.(4.3) is 



1 h 



{-l/2)X hh {a)g n °> 
(-!/2)(T 1 1 ^)- 



(4.2) 



[s/2] 



Vf(m,p) -> £>£(m,a,fc), 



(4.5) 



7>°(m,a,fc)= E ^Hm)pS;;;±..(a,fc). 

(i)GG 



play a role of boson fields and contractions, respectively. 



4.2 As a result, we come to the following general for- 
mulae. As far as homogeneous p-polynomials Vg{m,p) 
can be always represented as 

(i)€G 1 2 

i« G C, o= (4.3) 



In so doing, we are arrived to special j-degree ho- 
mogeneous polynomials in variables ( ^ ^ ) involved in 
primitive correlators, see Eq.(4.2). Polynomials of this 
type was introduced and named as hafnians by Caian- 
iello [66l [67] in the course of his QED investigations. 
Hafnians are the counterparts of phaffians and closely 
connected with permanents. The simplest nontrivial haf- 
nian ("1^°'^) of 2-degree is given below in the last 
three lines of Eq.(4.6). 



where the coefficients a s l (m) are functions in masses 
mi, I G £, it is sufficient to find the image of some gen- 
eral monomial entering into the sum over (i) G G in 
Eq.(4.3). Calculation according to the above mentioned 



4.3 Taking into account the very important applied 
meaning of an algorithm of constructing a family of 
homogeneous fc-polynomials "PfAm, a, k) from initial p- 
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polynomials Vf(m,p), we give some examples: 



Pr 

PlxPl'- 



n <Tl n CT2 „<T3 
Ph Ph Ph ~ 



p°Ip?A>u- 



= 0:1; 



lb 



j = l:H){X W2 ^ =:(££)}; 

i = i:(-i){(T 1 1 T 2 2 )[T|]+ 

+ V / x J 3 JL l 2 \ + I Z 2 *3 ^ h J J' 
*' /l ^2 ^3 ^4 * L ^1 ^2 ^3 ^4 -I ' 

j=i:(-l){(T 1 i r 2 2 )[T:r 4 4 ]+ 

+ Wi ; 3 ^ I i 2 k J + U, li ) ll 2 h J 

+ W 2 *3 I i Zl U J + U] 2 4 J ill «3 J 

+ (OTm 2 2 ]}' ( 4 - 6 ) 

J=2:(-i) 2 {(T^)(^) + 

, / fl °"3 \ / C2 CT4 \ / CTl <T 4 \ / <T 2 CT3 \ 1 

+ Ui ; 3 )\ h h ' + \ liU A Z 2 Z3 /J 

_. f lW <TlCT2 0'3Cr4 , i 

— ' V 2) Wiljlj h)- 



5. Parametric functions of FAs 

5.1 Now, let us formulate an algorithm of finding the 
parametric functions 

A(a), A(a,k), Yi(a,k), X w (a), 1,1' G £, 

of Feynman amplitudes. Of course, it is to be men- 
tioned that we can in principle use any one of the 
available approaches. Contributions to this subject 
have been made by many authors. We point out 
here the very incomplete list of quotes, namely pa- 
pers by Chisholm [55], Nambu [55], Symanzik [70] . 
Nakanishi [TJJ, Schimamoto [72], Bjorken and Wu [73], 
Peres [7J], Lam and Lebrun [75], Stepanov [75], Liu and 
Chow [77] , Cvitanovic and Kinoshita 78j , and books by 
Todorov [79], Speer [80], Nakanishi [ST], Zavialov [82], 
Smirnov [83) . in which many other citations can be 
found. Nevertheless, our algorithm seems to be very 
simple, but universal enough. It is named by the au- 
thor [84 , 85 , 32 as circuit-path algorithm. 

5.2 Suppose we have connected graph G(V, £ U £) 
with sets of vertices, V, of internal lines, £, of external 
lines, £, and with a certain relation of incidence between 
V and A = £ U £, described by an oriented incidence 



matrix en = [e VA ]u = 0, ±1, w, 6 V, I G A. In partic- 
ular, eu — 0, if the line I is nonincident to the vertex 
Vi] eu = 1, if the line I is outgoing from the vertex vf, 
eu = — 1, if the line / is incoming to the vertex tij. The 
fact that the set of all lines A is separated from the begin- 
ning into two mutually disjoint subsets £ and £ (incident 
properties of which are different) is very important both 
from algorithmic point of view and from potential pos- 
sibilities. In so doing, we need not replace here the set 
of external lines (incident to some one vertex) by some 
effective line, or assign the same orientation to all ex- 
ternal lines, as is usually done. Therefore, we can set 
the task of constructing parametric functions of whole 
graph via parametric functions of its subgraphs. As a 
result, circuit-path approach is naturally arose and the 
recursive structure of parametric functions of FAs has 
been obtained [55J |S5] . 

5.3 The set of external lines, £ , induces the single- 
valued decomposition of the set of all vertices, V, into 
the subset of external vertices, V ea:t , and the subset of in- 
ternal vertices, y m *. The set of internal lines, £, can be 
always decompose (as a rule in more than one way) into 
two mutually disjoint subsets, Ai and TV, which deter- 
mine some skeleton tree and corresponding to it co-tree 
subgraphs of the graph G. So, we have the following 
decomposition of the set A = £ U TV U A4 of all lines of 
the graph G into mutually disjoint subsets, £ , TV, and 
Ai. Then circuit-path algorithm requires the following 
steps: 

• Let us choose a subset AT C £ such that the sub- 
graph G(V, Ai U £ ), where Ai := C/Af, is a skeleton tree 
type graph and the subgraph G(V,Af U £) is a co-tree 
type graph. It is clear that this choice is ambiguous. It 
is shown in |84) . however, that the parametric functions 
are independent of any choice of Af. 

• Let us choose a vertex Vj S V which will be refered 
to as a basis vertex (or reference vertex, or zero point). 
It is clear that this choice is also ambiguous. But it is 
shown in [84] , that the parametric functions are again in- 
dependent of any given choice of Vj . From the viewpoint 
of practical calculations it seems reasonable to choose 
the basis vertex Vj as such a vertex to which the largest 
number of external lines of the graph are incident. 

• The choice of TV C £ and the basis vertex vj 
uniquely defines notions of basis circuits C(n), n G TV, 
and basis paths P(j\e), e G £ , namely. 

The basis circuit C{n) generated by the line n G Af 
is a union of the line n with the subset Ai(n) C Ai 
which forms a chain in Ai between vertices incident to 
the line n, i.e. C(n) := {n} U A4(n). The orientation in 
the circuit C(n) is defined by the orientation of the line 
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n g TV. 

The basis path P(j\e) generated by the line e G £ and 
the basis vertex Vj is a union of the line e with the sub- 
set A4(j\e) C M. which forms a chain in M. between a 
vertex incident to the line e G £ and the basis vertex vj , 
i.e. P(j\e) := {e} U A4(j\e). The orientation in the path 
P(j|e) is defined by the orientation of the line e G £. 

• By analogy with the incidence matrix e V A, that 
is more precisely can be refered to as the vertex- 
line incidence matrix, one introduces topologically the 
line- circuit e AA r [77J [75] HU [EH US], and the line-path 
EaeO) [SHI!!] incidence matrices, namely: 



[Caaa] In '■ — 

[e A£ (i)]i e : 



0, / g C(n), 
±1, Z G C(n); 

f 0, J£P(j|e), 
I ±1, /GP(j|e). 



(5.1) 



Here the plus or minus sign depends on whether the ori- 
entation of the line I G A coincide or not with the orien- 
tation of the circuit C(n) for e AjV - or the path P{j\e) for 
e A£ (j). As a result, the column-vector p A of all momenta 
pi, I G A, and submatrices of e AAf and e A£ (j), the rows of 
which are associated with the partition A = £ UTVU M, 
can be represented as follows [HH [85] : 



Pa =P A ext +ff *» VT = e A£ (j)k £ , pl nt = e KN p M - 
ees(j) = 1 £ £, e^ £ (j|AT) = Ojve, 
;(i -V) = 1 



= o. 



-£Af — u £AAj 



(5.2) 



A' — e {v/j}A4 e {v/j}Af- 



From now on, fc £ and p M are the column-vectors of 
external momenta fc e ,e G £, and independent inte- 
gration momenta p ni n G TV, respectively; AB is the 
| .4 1 x |,B | -rectangular matrix of zeros, and 1 AA is the 
\A 1 -dimensional unit matrix. Matrices e^ v/j}£ , e^ v/j}u , 



are submatrices of e VA . Their rows are defined 



by the set (V/vj) C V, and their columns are defined 
by the subsets £ , TV, M., respectively. The (|V| — 1)- 
dimensional square matrix e^ v/j)M is nonsingular, and 
detfe 

{v/]}m] — il- I n submatrices of the second and 
third lines of Eqs.(5.2), the subset TV is pointed out ex- 
plicitly, because of e M i e (j\N) ^ CV £ , and e M > £ (J\N) ^ 
e M e{jW) if TV' ^ TV, C = TV U M = TV' U M' , but 
e ££ (j\N) = e-eeijW') = X ££- 

• There exist the following very important "orthog- 



onality" relations [87l [84l [85]: 

e {v/j}A e AM = e {v/j}£. e £M — 0{v/j}AO 

[e V Ae A £(j)]ie =%[e(v*) £ ] e , 

e {v/i}A6A£0') = 0{v/j}£, 



(5.3) 



where e(v*) £ is the vertex- line incidence matrix of the 
"star" type graph G* :=< V*,£ > with the one vertex 
V* and the set of external lines £ of the graph G. The 
graph G* .—< V* , £ > is a result of shrinking of all ver- 
tices Vi € V, and all internal lines I G L, of the graph G 
to the single "black-hole" vertex V*. 

• On assigning to every internal line I G C the param- 
eter ai > 0, we define the circuit C A r J ^{a), path E ££ {j\a), 
and path- circuit Tl £A r(j\a) matrices |84 [ I85 | . according to: 



[G_^j^(a)] nn i : — [e^^a cc e cM\nn' — 

= ± J2 ai > 

leC(n)nC(n') 

[E ££ (j\a)] ee > := [el e (j)a cc e C£ (j)] ee/ 

leP(j\e)nP(j\e') 

[n £Ar (j|a)] e „:= [e£ £ (j)c\! CjC e £ jv]en = 
= ± J! ai - 

lGP(j\e)nC(n) 



(5.4) 



Here the plus or minus sign depends on the mutual ori- 
entations of the sets, over which the summation is per- 
formed, on their intersection. The plus sign corresponds 
to the case of coinciding orientations. It is clear that 
the explicit form of these matrices in any given case 
can be easily obtained by inspecting the graph. From 
now on a LL is the diagonal |£|-dimensional matrix, i.e., 

[otcc]ll> = Oil S W . 

• The parametric functions are derived by means of 
use the following matrices [BH [HI] : 

A ££ (j\a):= E ££ (j\a) - U £M (j\a)C^(a)U^(j\a), 
Y C£ (j\a):= e C£ (j) - e CJ ^C^(a}Il^(j\a), 
X cc (a):= e z;A rC A ; Af (a)eJ A ,, 

A(a):=detCW(a). (5.5) 

So, the quadratic A(a,k) and linear Yj(a, k), I G C, 
Kirchhoff forms in external momenta k e , e G £ , and the 
line-correlator functions Xu'(a), 1,1' G £, are defined 
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as [81186]: 

A(a,k):= (fcj • A e£ (j\a)k £ ) = 

= T,e t e'ee[ A ^(3\ a )]ee'{ke ' fee'), 

Yi(a,k):= Y ie (j\a)k £ = J^eesi^Ul 01 ^^, 

Y l£ {j\a) = e l£ (j) - e w C^(a)n^(j|a), 

X w {a) = e lM C^{a)e^ H . (5.6) 

Here e !jV and e l£ (j) are row-vectors of matrices (5.1)- 
(5.2) for the line I G C. 

5.4 It should be mentioned that functions A (a) and 
A{a, k) do not depend on the orientation of internal 
lines. However, when the orientation of the line / is 
changed, parametric functions Yi(a,k) and Xui(a) are 
reversed their sign. 

It is also useful to represent quantities A(a, k) and 
Y £ {a, k) in a form exhibiting a special role of matrices 
X cc (a) and^(a) [TEUM]: 

A(a, k) = 

= {pt xt {k) T -[a cc - a cc X cc {a)a cc ]pl xt {k)), 
p?\k) = e C£ (j)k £ , p e £ xt (k) = k £ , p e f(k)=0^, 
Y c (a,k) = [l cc - X cc {a)a cc ]pl xt {k) = 
= p™ t (k)-Yj nt (a,k), 
Y £ mt (a,k):=X ££ (a)a ££ p e £ xt (k) 7 
X ££ (a):= e £M X NM (a)e £H: 

X UM {a):=C^{a), (5.7) 

where k £ is the column- vector of external momenta fc e , 
e G £ . There hold the following homogeneous properties: 

A(pa) =p lcl A(a), X w (pa) = p- l X w (a), 
A{pa 1 rk) — pr 2 A(a, k), Yi(pa, rk) — rYi(a, fc), 
7>°(m, pa, rk) = P - 3 t s - 23 P°-(m, a, fc). (5.8) 

5.5 Now we exhibit some important properties of a- 
parametric functions |45j . Let us introduce quantities, 

K ££ -=X ££ a ££ , K cc := a ££ X ££ , 

V ££ :=\ ££ -K\ £1 i = r,l. (5.9) 

Using Eqs.(5.3)-(5.5), we find that matrices K ££ (a) and 
L ££ [a) are projectors with next properties: 



- 11 r.r.- 1 ^ cc 



K " cci L ££ L ££ — L ££ , i — r,l, 



K l £.cL ££ — ££ , K cc a ££ L r ££ — ££ , 



(5.10) 



From Eqs.(5.10), we get some relations between products 
of X cc , a ££ , and Y C£ : 

(X ££ a ££ j X ££ = X ££ a ££ X ££ — X cc , 

iPcc) X ££ — CC , 

(X ££ a ££ ) Y ££ = X ££ a ££ Y ££ = ££ , (5-11) 
(L^) m Y ££ = Y ££l 

Tr[(Ki c ) m ] = Tr[K cc ) = £ ie£ a,X„(a) = M> 

Tr[(i'U m ] =1V[LL] =\M\, i = r,l, 

and the following relations between quadratic A(a, fc) 
and linear Yi(a, fc), I G C, Kirchhoff forms: 

A(a, fc) = (Y C T ■ a ££ p ext ) = {p extT ■ a ££ Y £ ) = 

^ZiecMPi ext (k)-Y l (a,k))= (5.12) 
= (Y £ T ■ a ££ Y £ ) S £ ie£ a z lf(a,fc). 

There exist also the following relations, 

e v£ k £ + e v£ Y £ (a, fc) = V , e £ ^a ££ Y £ (a, fc) = 0^, 
e v£ X ££ (a) = V£ , 

K lc e ce(j) = -ecxY M£ {j\a) = e £M K r M£ e ££ (j), 
K r ££ e cN = e c.Mi (5.13) 
{Y c mtT ■ a ££ Y £ ) = 

= (Y £ mtT ■ a ££ pl xt ) - (Y £ mtT ■ a ££ Y £ mt ) = 0. 

Two relations in the first line of Eqs.(5.13) in our case of 
a-parametric functions are analogs of the first and the 
second Kirchhoff laws in electric networks. Similarly, in 
the third line of Eqs.(5.12) we find in our case an analog 
of the well-known expression for a power dissipated in 
electric networks. 

The present author wishes to express his gratitude 
thanks to Reviewer for constructive review and valuable 
suggestions that enhances the improving of the article. 
The paper is based on the report presented at the Bo- 
golyubov Kyiv conference "Modern Problems of The- 
oretical and Mathematical Physics", September 15-18, 
2009, Kyiv, Ukraine. 
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